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A  MONTE  CARLO  MODEL  FOR  DETERMINING  COPPERHEAD  PROBABILITY 
OF  ACQUISITION  AND  MANEUVER 


1.  INTRODUCTION 

AMSAA's  Probability  of  Acquisition  and  Maneuver  (PAM)  model  has 
proven  to  be  a  useful  tool  for  evaluating  the  COPPERHEAD  weapon  system. 
Different  versions  of  the  model  have  been  used  in  two  different  ways. 

A  "stand  alone"  version  of  the  model  was  used  in  the  COPPERHEAD 
analysis  documented  in  Reference  1.  This  analysis  evaluated  the  sensitivity 
of  COPPERHEAD  system  performance  to  a  large  number  of  factors,  including: 

t  cloud  ceiling 

t  designator  power 

•  designator  range 

t  target  reflectivity 
t  target  location  error 

•  atmospheric  transmission 
t  seeker  sensitivity 

t  gun-target  range 

•  unguided  delivery  error 

An  improved  version  of  the  stand  alone  PAM  model  was  later  used 
to  evaluate  the  sensitivity  of  COPPERHEAD  system  performance  to: 

•  designator-target-howitzer  azimuth  angle  (ANGLE-T) 
t  target  velocity 

t  response  time 
t  additional  delay  time 
t  target  heading  angle 

•  point  of  target's  closest  approach  to  Predicted  Intercept 
Point  (PIP) 


IChernick,  Julian  A.,  Richard  C.  Scungio,  Michael  Starks,  Utility  of 
COPPERHEAD  With  Ground  Laser  Designation  in  a  European  Battlefield 

Environment  (U),  Technical  Report  No.  257,  US  Army  Materiel  Systems 

Analysis  Activity,  Aberdeen  Proving  Ground,  MD,  December  1978, 
(CONFIDENTIAL  report). 


5 


The  results  of  this  analysis  were  published  as  part  of  the  COPPERHEAD 
COEA  (Reference  2).  Other  "stand  alone"  uses  to  which  modified  versions 
of  PAM  have  been  applied  include  data  generation  for  the  Advanced  Anti- 
Armor  Vehicle  Evaluation  (ARMVAL)  tests,  and  for  a  forthcoming  Joint 
Technical  Coordinating  Group  for  Munitions  Effectiveness  (JTCG/ME)  publi¬ 
cation  on  COPPERHEAD. 

In  addition  to  its  use  as  a  stand-alone  model,  PAM  serves  as 
one  of  the  preprocessors  for  AMSAA's  COPPERHEAD  Operational  Performance 
Evaluation  (COPE)  model.  For  a  description  of  the  COPE  model,  including 
details  of  that  model's  interface  with  PAM,  see  Reference  3. 


Thus  far,  the  model  applications  mentioned  have  concerned  the 
COPPERHEAD  system.  First-order  performance  estimates  for  other  weapon 
systems  have  also  been  developed  through  use  of  modified  versions  of  the 
PAM/COPE  models.  These  systems  include  HELLFIRE  and  extended  range 
COPPERHEAD  (Reference  4). 


The  purpose  of  this  report  is  to  document  the  structure  of  the 
PAM  model  so  that  other  activities  may  more  easily  use  it  in  related 
analyses.  The  report  is  organized  as  follows: 

•  The  assumptions  made  in  the  course  of  constructing  the 
model  are  discussed. 


•  A  general  overview  of  the  model  structure  is  presented. 

•  The  way  in  which  the  acquisition  portion  of  the  COPPERHEAD 
trajectory  is  modeled  is  described  in  detail. 

•  The  way  in  which  the  maneuver  portion  of  the  COPPERHEAD 
trajectory  is  modeled  is  described  in  detail. 

•  Appendix  A  lists  the  inputs  required  to  drive  the  model, 
along  with  the  appropriate  units. 

•  Appendix  B  contains  a  copy  of  the  FORTRAN  SOURCE  LIST. 

•  Appendix  C  gives  a  sample  case  with  input  and  output. 


^Cost  and  Operational  Effectiveness  Analysis  (COpPERHEAD,  COEA) (U ) , 

ACN  16812,  US  Army  Field  Artillery  School,  FT  Sill,  OK,  October  1979, 
(SECRET  report). 

-^Sandmeyer,  Richard  S. ,  COPE  Computer  Program:  User  and  Analyst  Manuals, 
US  Army  Materiel  Systems  Analysis  Activity,  Aberdeen  Proving  Ground,  MD, 
Technical  Report  to  be  published. 

^Chernick,  Julian  A.,  Preliminary  Analysis  of  Extended  Range  COPPERHEAD 
Operational  Performance  (U),  GWD  Interim  Note  No.  G-85,  US  Army  Materiel 
Systems  Analysis  Activity,  Aberdeen  Proving  Ground,  MD,  January  1980, 
(CONFIDENTIAL  report). 
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2.  ASSUMPTIONS  AND  LIMITATIONS 


As  the  name  of  the  model  suggests,  the  PAM  model  computes  the 
Probability  of  a  COPPERHEAD  projectile  being  able  to  both  Acquire  a  target 
Fy  sensing  reflected  laser  energy,  and  Maneuver  to  that  target  once  it 
has  been  acquired.  The  fundamental  difference  between  the  probability  of 
acquisition  and  maneuver  (PAM)  and  probability  of  hit  (Pj-j) ,  is  that 
estimates  of  Pj-i  for  COPPERHEAD-type  systems  must  include  the  effect  of 
laser  energy  overspill  and  underspill. 

Such  effects  are  a  complex  function  of  the  entire  time-history 
of  laser  pulses,  and  simulating  those  effects  requires  detailed  modeling 
of  the  system's  seeker  logic.  The  Laser  Designator  Weapon  System  Simula¬ 
tion  (LDWSS)  model  (Reference  5)  does  simulate  the  seeker  logic,  so  the 
resulting  estimates  of  P^  include  the  effects  of  laser  energy  overspill 
and  underspill. 

The  PAM  model  does  not  simulate  these  effects;  however,  there 
is  reason  to  believe  that  this  limitation  is  not  too  severe.  Section  2.3.2 
of  Reference  6  presents  LDWSS  estimates  of  Ph  for  the  COPPERHEAD  system 
under  various  conditions.  Under  conditions  of  high  visibility,  high  cloud 
ceiling,  low  errors,  and  GLLD  designator,  it  is  plausible  to  assume  that 
any  degradation  in  P^|  is  due  to  spill  over/spill  under.  As  the  data  shows, 
there  is  little  degradation  in  P^  against  a  fully  exposed  moving  target 
out  to  3  km.  Therefore  there  is  little  problem  with  spil lover/spil 1  under 
against  such  a  target  out  to  3  km.  However,  for  a  partially  exposed 
target  or  a  target  at  longer  range,  the  spil  lover/spil  under  problem  is 
more  severe.  Under  these  conditions  the  probability  of  acquisition  and 
maneuver  is  a  poor  estimator  of  probabil ity  of  hit. 

A  second  limitation  of  the  PAM  model  is  the  use  of  a  lambertian 
reflectance  distribution  (cosine  law)  of  energy  from  the  target  rather 
than  specific  reflectivity  maps  generated  from  a  three-dimensional  target 
description.  While  significant  differences  could  exist  in  terms  of  the 
actual  shape  of  the  acquisition  volume  for  each  reflected  laser  pulse, 
the  spot  jitter  and  the  time-variabil ity  of  target  heading  is  probably 
sufficient  to  smooth  out  the  shape  of  the  acquisition  volume  in  such  a 
way  that  the  cosine  law  is  approximately  correct. 

Additional  limitations  arise  because  the  model  uses  Monte  Carlo 
sampling.  The  results  are  somewhat  noisy  (0.0  percent)  when  a  sample  size 
of  100  is  used;  trends  which  intuitively  should  be  monotonic  are  not 
always  so.  Moreover,  when  larger  sample  sizes  are  used,  the  model  becomes 
fairly  time  consuming  to  run.  Still,  results  of  good  quality  can  be 
obtained  with  greater  ease  and  at  lesser  expense  than  by  use  of  LDWSS; 
the  LDWSS  model  also  uses  Monte  Carlo  sampling. 


5Lewis,  C.L.,  A.G.  Nichols,  and  A.W.  Lee,  User's  Guide  for  the  Phase  I 
Laser  Designator/Weapon  System  Simulation  (LDWSS)  of  the  COPPERHEAD 

Guided  Projectile  System, Vol  I,  Technical  Report  RG-/7-2S.  US  Army 

Missile  Command,  Redstone  Arsenal,  AL,  July  1977,  (UNCLASSIFIED  report). 
^Independent  Evaluation  Report  for  the  155mm  XM712  COPPERHEAD  (U),  IER 
No.  6-80,  US  Army  Materiel  Systems  Analysis  Activity,  Aberdeen  Proving 
Ground,  MD,  July  1979,  (CONFIDENTIAL  report). 
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3.  GENERAL  OVERVIEW 


The  overall  flow  of  the  main  PAM  model  is  shown  in  Figure  1. 

The  most  important  part  of  the  model  is  devoted  to  answering  the  two 
questions  which  appear  near  the  center  of  Figure  1:  "Does  the  projectile 
acquire?"  and  if  so,  "Can  the  projectile  maneuver  to  the  target?" 

Figure  2  shows  the  general,  situation  of  a  ground-based  laser 
designating  a  target.  The  acquisition  volume  is  the  volume  within  which 
there  is  sufficient  reflected  laser  energy  for  the  COPPERHEAD  projectile 
to  acquire  the  target.  The  optical  properties  of  the  atmosphere  determine 
the  extent  of  laser  energy  attenuation  along  both  the  laser-target  path 
and  the  target  projectile  path.  The  cloud  ceiling  acts  as  an  energy 
cutoff  level,  prohibiting  acquisition  until  the  projectile  descends 
below  the  cloud  layer.  Cloud  cover  is  treated  as  opaque;  beneath  a 
cloud  layer  the  visibility  is  assumed  uniform. 

Once  the  COPPERHEAD  projectile  breaks  through  a  cloud  ceiling, 
it  acquires  a  target  only  if  reflected  laser  energy  reaches  the  seeker  in 
sufficient  quantity.  This  process  is  illustrated  in  Figure  3. 

Since  there  is  a  ceiling  at  altitude  A,  acquisition  is  impossible 
above  this  ceiling.  At  altitude  B,  acquisition  is  possible  but  does  not 
take  place  because  insufficient  energy  reaches  the  seeker.  At  altitude  C, 
sufficient  energy  reaches  the  seeker  for  acquisition  to  take  place. 

The  (x,y)  coordinates  of  intersection  with  the  acquisition 
volume  at  a  given  altitude  plane  is  a  function  of  unguided  delivery  error. 
Given  (x,y,z)  coordinates  of  acquisition,  the  limits  of  projectile  maneuver 
in  the  ground  plane  are  fixed.  For  the  case  illustrated  in  Figure  3,  the 
target  is  within  the  limits  of  projectile  maneuver,  so  the  engagement  is 
a  success. 
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FIGURE  1  CONCEPTUAL  FLOWCHART  OF  PAM  MODEL 
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Figure  2  Laser  Designating  a  Target. 
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Copperhead  Acquisition  and  Maneuver. 


4.  ACQUISITION  METHODOLOGY 


The  equations  used  in  PAM  to  describe  the  laser  energy  transmis¬ 
sion  are  the  same  as  those  used  in  the  LDWSS  model  (Reference  5).  The 
laser  beam  energy  signal -to-threshold  (S/T)  ratio  is: 


S/T  = 


EdTdTspcos0 


where 

E.J-  =  threshold  energy  density  at  the  seeker  aperature  (J/km^) 

E,j  =  laser  designator  energy  (J) 

T,j  =  designator-to-target  transmission  coefficient 
Ts  =  target-to-seeker  transmission  coefficient 
p  =  target  reflectivity 

■  0  =  lambertian  angle  (angle  from  seeker  LOS  to  designator  beam) 

Rs  =  slant  range  from  target  to  seeker  (km) 

The  transmission  coefficients  of  the  laser  equation  are  complex  functions 
of  altitude,  general  atmospheric  condition,  and  wavelength.  These  coef¬ 
ficients  are  calculated  as  functions  of  visibility,  altitude  (Hs), 
projectile  range  to  the  target  (Rs)  and  designator  range  to  the  target 
(R<j)  in  km:* 


Td  =  e"  Rd 

.  _  e-.00025H$ 

e  "T  .00025H 
Ts  =  e  -YRS 


Rs’  for  Hs  >  0 
,  for  Hs  <  0  (Ft) 


*It  can  be  seen  from  examination  of  these  two  expressions  that  PAM 
and  LDWSS  both  assume  that  electro-optical  transmissivity  improves  as 
a  function  of  increasing  altitude.  Recent  work  done  at  the  Atmospheric 
Sciences  Laboratory  indicates  that  this  assumption  is  not  always  true 
(Reference  7).  During  certain  conditions  of  fog  and  haze,  transmissivity 
may  be  as  much  as  two  orders  of  magnitude  worse  at  150m  above  the  ground 
than  at  ground  1 evel . 

'Pinnick,  R.G.,  et.al..  Vertical  Structure  in  Atmospheric  Fog  and  Haze 
and  Its  Effects  on  IR  Extinction,  Atmospheric  Sciences  Laboratory, 

White  Sands  Missile  Range,  NM,  EC0M-TR-0010,  July  1978. 
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The  atmospheric  attenuation  coefficient  (y)  is  a  function  of 
vi sibil i ty  (VIS) : 

v  _  .0019(.519)Q 
1  VIS 


with  Q  determined  as  a  function  of  visibility  in  kilometers; 
1/3 

VIS  ,  0<VIS<6 
5 


VIS 

0.86  +  3q-  ,  6_<VIS<9 


VIS 

0.98  +  5Q-  ,  <VIS<12 


VIS 

!.!5  +  200  »  VIS>12 


The  resultant  visibility  volume  around  the  target  has  its 
maximum  length  along  the  direction  of  the  designator  and  is  of  negligible 
extent  for  angles  greater  than  90  degrees  from  the  designator-target  line. 

The  model  uses  a  nominal  input  value  for  angle  T  (FO-Target- 
Howitzer  azimuth  angle).  However,  the  cosine  law  of  reflectance  is  not 
applied  directly  to  that  angle,  but  to  the  input  angle  adjusted  for  the 
actual  target  location  present  in  a  particular  Monte  Carlo  replication 
after  unguided  errors  and  target  location  errors  are  sampled. 
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5.  MANEUVER  METHODOLOGY 


If,  for  a  given  Monte  Carlo  replication,  the  model  determines 
that  a  COPPERHEAD  projectile  acquires  a  target,  then  computations  are  made 
to  determine  whether  the  projectile  can  maneuver  to  that  target.  This  is 
accomplished  by  means  of  maneuverability  footprints.* 

For  a  given  gun-target  range,  mode  of  fire,  angle  of  fall,  and 
altitude  at  which  initial  acquisition  takes  place,  the  footprints  circum¬ 
scribe  an  area  in  the  ground  plane  within  which  a  reliable  projectile 
can  successfully  maneuver.  This  area  is  the  intersection  of  the  seeker 
f i el d-of-view  projected  into  the  ground  plane  and  the  extreme  limits  of 
projectile  maneuver  capability.  Because  the  footprints  lack  radial 
symmetry,  they  are  input  as  a  series  of  distances  as  a  function  of  angle 
from  the  predicted  target  intercept  point  (PIP). 

The  model  considers  three  kinds  of  error  source:  unguided 
delivery  errors,  random  target  location  errors,  and  bias  target  location 
errors.  Given  values  for  these  three  error  terms,  the  model  determines 
whether  the  target  is  in  the  footprint  at  the  time  of  round  arrival. 

Unguided  delivery  errors  are  associated  with  the  part  of  the 
COPPERHEAD  trajectory  between  launch  and  acquisition.  Such  errors  have 
the  effect  of  shifting  the  location  of  the  footprint  in  the  ground  plane; 
they  are  monte  carlo  sampled  for  each  simulated  trajectory. 

Standard  deviations  for  these  errors  in  range  and  deflection 
are  shown  in  Figure  4.  The  information  shown  was  supplied  by  Martin 
Marietta  Corporation,  and  was  generated  by  use  of  six-degree-of-freedom 
simulation  techniques.  Range  errors  are  larger  than  deflection  errors 
primarily  because  of  COPPERHEAD'S  relatively  shallow  angle  of  fall. 
Although  the  range  and  deflection  errors  are  not  too  different  in  the 
plane  normal  to  the  velocity  vector,  when  they  are  projected  into  the 
ground  plane,  the  range  error  becomes  elongated. 

*Both  ARRADCOM  and  Martin  Marietta  Corporation  have  supplied  AMSAA  with 
these  footprints.  For  details  on  the  ARRADCOM  model  which  generates 
footprints  see  Reference  8. 

Amoruso,  Michael,  J.,  Tice  F.  DeYoung,  Dennis  D.  Ladd,  and  Roger  D. 

Schulz,  A  Comprehensive  Digital  Flight  Simulation  of  the  Cannon 
Launched  Guided  Projectile,  Rodman  Laboratory,  Rock  Island,  IL, 

January  1977,  R-TR-77-007  (UNCLASSIFIED  report). 
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Figure  4.  Unguided  Delivery  Errors 


The  random  component  of  target  location  error  (a)  can  be  input 
to  PAM  directly  or  computed  internally  according  to  the  following  algorithm 
derived  from  unpublished  work  of  Julian  Chernick.*  The  algorithm 
contains  two  error  terms  and  a  parameter: 

c2TLE  =  a2 FO  +  (a2S  x  T2) 

where 


a  TLE  =  variance  of  target  location  error 
2 

o  FO  =  variance  of  FO/FDC  error  in  locating  FO  position 
2 

o  S  =  variance  of  FO  error  in  estimating  target  speed 
T  =  total  system  response  time 

The  time  term  in  the  algorithm  is  the  sum  of  the  expected  re¬ 
sponse  time  and  the  unanticipated  delay  time.  Expected  response  time  is 
the  average  time  between  the  beginning  of  the  FO's  call  for  fire  and  the 
time  of  round  arrival  on  target;  it  is  an  input  to  the  model.  Unantici¬ 
pated  delays  are  played  parametrically  in  the  model  with  values  of  0,  30, 
90,  150  and  300  seconds. 

The  random  TLE  algorithm  above  allows  aTLE  to  be  computed 
about  a  single  target  vehicle,  for  either  a  preplanned  target  or  a  target 
of  opportunity.  Since  the  PAM  model  was  designed  to  evaluate  COPPERHEAD 
against  groups  of  target  vehicles  as  well  as  against  a  single  target 
vehicle,  a  method  was  derived  to  generate  oTLE  to  the  nearest  target 
vehicle  when  more  than  one  vehicle  is  present  in  a  target. 

Based  on  work  reported  in  Reference  9,  the  following  relation¬ 
ship  is  used  for  computing  oTLE  when  the  PIP  is  bracketed  by  target 
vehicles : 


oTLE  (multiple)  =  .68  oTLE  (single) 


*A  more  general  account  of  this  type  of  error  has  recently  been  developed 
by  Chernick  (Reference  10).  The  numerical  results,  however,  are  similar 
to  those  resulting  from  the  present  algorithm. 

9Weaver,  Jonathan  M. ,  and  Lawrence  Bowman,  Multiple  Target  Simulation 
(MUTSI)  -  A  Discrete  Monte  Carlo  Technique  That  Evaluates  the  Availability 

of  Multiple  Enemy  Ground  Targets  in  a  GLLD/COPPERHEAD  Target  of  Opportunity 

Situation,  GWD  Interim  Note  No.  G-61 ,  US  Army  Materiel  Systems  Analysis 

Activity,  Aberdeen  Proving  Ground,  MD,  July  1979  (UNCLASSIFIED  report). 
lOChernick,  Julian  A.,  Moving  Target  Location  Errors  for  Ground  Targets, 

US  Army  Materiel  Systems  Analysis  Activity,  Aberdeen  Proving  Ground,  MD, 
September  1980,  to  be  published. 
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The  final  error  source  considered  in  "he  PAM  model  is  a  bias 
TLE.  The  random  TLE  discussed  above  is  computed  on  the  assumption  that 
the  mean  target  position  at  round  impact  is  the  PIP.  This  assumption 
will  seldom  be  satisfied  in  combat,  hence  the  resulting  offsets  from  the 
PIP  are  treated  as  bias  TLEs. 

One  important  source  of  bias  error  is  unanticipated  delays  in 
the  time  required  to  get  a  COPPERHEAD  round  on  -arget.  If  the  FO  estimates 
it  will  take  100  seconds  to  get  a  round  on  target,  but  it  actually  takes 
200  seconds,  then  the  target  may  overrun  the  footprint  before  the  round 
arri ves. 


In  addition  to  the  contribution  of  tine  delays  to  bias  TLE, 
there  are  factors  which  could  cause  the  target's  point  of  closest  approach 
(PCA)  to  the  PIP  to  differ,  on  the  average,  fron  zero.  If  the  footprint 
(aimpoint)  is  preplanned,  it  would  be  unreasonable  to  expect  that  potential 
targets  would  be  headed  directly  towards  the  PIP.  And  for  a  target-of- 
opportunity,  there  is  a  possibility  of  large  changes  of  direction  after 
the  command  to  fire  is  given. 

Bias  error  is  played  in  PAM  as  follows:  If  a  COPPERHEAD 
projectile  approaches  the  target  from  the  direction  of  negative  y,  the  x 
and  y  components  of  bias  TLE  are  given  by 

XBIAS  =  PCAX  +  S i n ( H ) V  x 

YBIAS  =  PCAy  -  Cos(H)Vt 


where 


PCA  -  point  of  closest  approach 
H  -  target  heading  angle 
V  -  target  velocity 
t  -  unanticipated  delay 

After  these  computations  are  made,  the  random  and  bias  TLEs  are 
summed  to  yield  the  target's  location  on  the  ground.  The  distance  from 
that  location  to  the  PIP  is  computed,  and  is  conpared  to  the  distance 
from  the  PIP  to  the  edge  of  the  maneuverability  footprint  for  an  equivalent 
angle.  If  the  target  is  within  the  footprint,  then  the  replication  counts 
as  a  success;  otherwise  it  does  not. 

6.  SUMMARY 

This  report  presents  the  general  structure  of  the  PAM  model,  and 
describes  the  modeling  of  the  acquisition  and  maneuver  positions  of  the 
COPPERHEAD  trajectory  in  detail.  In  addition,  a  description  of  the  input 
variables,  a  FORTRAN  source  listing,  and  a  sample  case  are  presented. 
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APPENDIX  A 


PAM  INPUTS 

Table  A-l  provides  a  list  of  the  inputs  required  for  the  PAM 
model.  These  are  read  in  with  list-directed  read  statements  as  indicated 
in  the  FORTRAN  listing  in  Appendix  B. 

In  addition  to  the  input  variables  which  are  read  in,  three 
arrays  are  filled  in  the  program  and  the  values  could  be  changed  at  the 
discretion  of  the  program  user.  These  arrays  are: 

GAMARY  (11) 

H  (6) 

DTS  (5) 

The  GAMARY  array  contains  the  laser  energy  attenuation  coef¬ 
ficients,  as  a  function  of  atmospheric  visibility  in  km,  from  1  to  11 
km. 

The  H  array  defines  the  heights  at  which  the  check  for  acquisi¬ 
tion  is  made,  from  highest  to  lowest,  in  meters. 

The  DTS  array  holds  the  values  of  unanticipated  delay  times 
played  in  seconds,  as  discussed  in  the  Maneuver  Methodology  section  of 
this  report. 
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TABLE  A-l  PAM  INPUT  VARIABLES 


VARIABLE  NAME 

MEANING 

UNITS 

IMF 

Mode  of  f i re 

1  =  preplanned 

2  =  target  of 

opportunity 

ETH 

Seeker  Energy  Threshold 

Joules/ SQKM 

AOF 

Fly  under  Fly  out  (FUFO) 

Angle  of  Fall 

Degrees 

TH 

Target  Heading  Angle 

Degrees 

PCA 

Point  of  Closest  Approach 

Meters 

AZDT 

Nomi nal  Angl e  T 

Degrees 

V 

Target  Velocity 

M/S 

RHO 

Target  Reflectivity 

N/A 

ED 

Designator  Energy 

Joul es 

TR 

Nominal  Response  Time 

Including  TOF 

Seconds 

NK 

Number  of  Monte  Carlo  Cases 

N/A 

RNG 

Gun  Target  Range 

Meters 

ACCX 

Bal 1 i stic  Error  x  ( a) 

Meters 

ACCY 

Bal 1 i stic  Error  y  ( a) 

Meters 

IDRMN 

Minimum  Designation  Rng 

Km 

IDRMX 

Maximum  Designation  Rng 

Km 

I VMX 

Maximum  Visibility  Rng 

Km 

NI(J) 

#  of  Pts  in  which  Jth 

Footprint  is  input 

N/A 

THEMN(I,J) 

Ith  Angle  from  PIP  in  Jth 
footpri nt 

Degrees 

DISMH(I,J) 

Ith  Distance  from  PIP  to 
edge  of  Jth  Footprint 

Meters 
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APPENDIX  B 


FORTRAN  SOURCE  LIST 

This  Appendix  contains  a  FORTRAN  listing  of  the  PAM  model  as 
configured  for  a  CDC  7600  computer.  The  version  of  the  program  given  here 
is  for  interface  with  the  COPE  model;  comment  cards  at  the  beginning  of 
the  listing  indicate  necessary  deletions  for  use  as  a  stand-alone  model. 

Two  subprograms,  generally  available  as  system  routines  (for 
example,  on  the  Ballistic  Research  Laboratory  computer  at  APG)  are 
included  to  facilitate  program  portability.  These  are  subroutines  NRAN31 
which  generates  pairs  of  normally  distributed  random  numbers,  and  sub¬ 
routine  DVDINT,  which  does  divided  difference  interpolation. 
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APPENDIX  C 


SAMPLE  CASE 

A  list  of  input  values  as  used  in  a  sample  PAM  run  is  shown  in 
Table  C-l.  These  values  are  for  the  purpose  of  illustration  only.  The 
values  are  read  into  the  input  variables  in  the  order  that  those  variables 
are  listed  in  Table  A-l.  For  example,  the  first  card  indicates  that  the 
program  run  is  for  a  preplanned  footprint  situation;  the  second  card 
indicates  that  the  seeker  sensitivity  being  played  is  .00003  Joules/Km^-. 

Part  of  the  output  which  resulted  from  running  the  program  with 
the  sample  input  stream  is  shown  in  Figure  C-2.  For  each  combination  of 
designation  range  and  unanticipated  delay  time  there  is  a  6  x  10  output 
matrix.  In  addition,  for  the  case  where  there  is  no  unanticipated  delay, 
probabilities  are  printed  for  both  single  and  multiple  target  TLEs. 

Output  values  are  configured  in  each  matrix  as  follows: 

Highest  ceiling  .  Lowest  Ceiling 

Visibility  =  1  Km 


I 


Vi  si bi 1 ity  =  10  Km 


Variable 

(See  Table  A1 )  Sample  Values (s) 


IMF 

1 

ETH 

.000030 

AOF 

20.0 

TH,PCA 

0.0. 0.0 

A2DT,V,RH0,ED,TR 

25. 0.5.,. 10,. 1,106. 

NK,RNG,ACCX,ACCY 

100,8000. ,73. ,366. 

IDRMN,1DRMX,IVMX 

1,7,10 

NI(1) 

8 

THEMN(1 ,8) 

0.0. 20. ,30. ,45. ,60. ,90. 

,135. 

,180. 

DISMH(1 ,8) 

1200., 1  100., 1200., 1000. 

,600. 

,500. ,400. ,700 

NI(2) 

8 

THEMN(2,8) 

0.  ,20.  ,30. ,45. ,60.  ,90. , 

135., 

180. 

DISMH(2,8) 

1200.,  11 00., 1200., 1000. 

,600. 

,500.  ,400. ,700 

NI(3) 

8 

THEMN(3,8) 

0., 20., 30., 45., 60. ,90., 

135. , 

180. 

DISMH(2,8) 

1100. ,1000. ,1100. ,900., 

500. , 

400. ,300. ,700. 

NI  (4) 

8 

THEMN(4,8) 

0.  ,20.  ,30. ,45., 60. ,90. , 

135.. 

180. 

DISMH(4,8) 

1000. ,900. , 10C0. ,800. ,500. .300. .300. .500. 

NI  (5) 

8 

THEMN(5,8) 

0., 20., 30., 45., 60., 90., 

135., 

180. 

DISMH(5,8) 

600. ,500. ,600. ,400. ,500 

.,400 

.  ,400. ,600. 

NI  (6) 

8 

THEMN(6,8) 

0.  ,20.  ,30. ,45. ,60. ,90. , 

135., 

180. 

DISMN(6,8) 

400. ,300. ,400. ,300. ,400 

.,300 

.  ,300.  ,400. 

FIGURE  C2  SAMPLE  PAM  INPUTS 
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PROBABILITIES  OF  ACQUISITION  AND  MANEUVER 


Visibility  Designation  Range  1  km 

Range  Time  Delay  30s 

(km)  Cloud  Ceiling  (Ft) 


4500 

3000 

2500 

2000 

1500 

1000 

1 

.20 

.20 

.20 

.20 

.20 

.16 

2 

.48 

.48 

.46 

.44 

.44 

.24 

3 

.46 

.46 

.40 

.34 

.31 

.22 

4 

.54 

.54 

.51 

.49 

.47 

.30 

5 

.54 

.54 

.49 

.47 

.45 

.28 

6 

.54 

.54 

.49 

.44 

.39 

.25 

7 

.51 

.51 

.47 

.44 

.42 

.29 

8 

.54 

.54 

.50 

.48 

.48 

.28 

9 

.45 

.45 

.40 

.36 

.33 

.21 

10 

.55 

.55 

.48 

.45 

.42 

.25 

FIGURE  C2  SAMPLE  PAM  OUTPUT 
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